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Resume : Nous proposons une methode de reduction de modele, pour des ecoulements 
actionnes, basee sur la regularisation d'un probleme inverse. Le probleme inverse vise a 
minimiser l'erreur entre les predictions du modele et des simulations de reference. Les 
parametres a identifier sont les coefficients d'une expansion polynomiale qui modelise 
les dynamiques temporelles d'un nombre limite de modes globales. Ces modes sont 
obtenus par Decomposition Orthogonale aux valeurs Propres (POD). II s'agit d'une 
methode pour calculer les elements les plus representatifs, en termes d'energie, d'une 
base de donnees de simulations. II est montre que des modeles bases sur une simple pro- 
jection de Galerkin et sur des techniques classiques de calibration ne sont pas viables. 
lis sont soit mal poses, soit donnent une approximation erronee de la solution des qu'ils 
sont utilises pour predire des dynamiques n'appartenant pas a la base de donnees ori- 
ginale. En revanche, des evidences numeriques montrent que la methode que nous pro- 
posons permet de construire des modeles robustes, qui reagissent correctement a des 
variations de parametres, et pourraient done tre utilises pour des problemes de controle 
d'ecoulement. 

Mots-cles : modeles reduits, controle, problemes inverses 
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Robust model identification of actuated vortex wakes 



Abstract: We present a low-order modeling technique for actuated flows based on the 
regularization of an inverse problem. The inverse problem aims at minimizing the error 
between the model predictions and some reference simulations. The parameters to be 
identified are a subset of the coefficients of a polynomial expansion which models the 
temporal dynamics of a small number of global modes. These global modes are found 
by Proper Orthogonal Decomposition, which is a method to compute the most repre- 
sentative elements of an existing simulation database in terms of energy. It is shown that 
low-order control models based on a simple Galerkin projection and usual calibration 
techniques are not viable. They are either ill-posed or they give a poor approximation 
of the solution as soon as they are used to predict cases not belonging to the original 
solution database. In contrast, numerical evidence shows that the method we propose 
is robust with respect to variations of the control laws applied, thus allowing the actual 
use of such models for control. 

Key-words: reduced order models, control, inverse problems 
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1 Introduction 

We consider the problem of describing the dynamics of an infinite dimensional 
system using a small number of degrees of freedom. In particular we concentrate on 
the problem of devising accurate and robust models of actuated fluid flows past bluff 
obstacles. These flows are dominated by the presence of large-scale vortices due to 
massive separation, and are good candidates for a low-dimensional representation. The 
point of view that we privilege is empiric : the functional space in which we seek 
the low-dimensional solution is derived using proper orthogonal decomposition (POD) 
lfl6ll . POD makes use of simulation databases to determine optimal functional spaces 
in terms of solution representation. A vast literature concerning this way of modeling 
fluid flows exists JBBSB^ and some results show the possible interest of using POD 
in applications such as flow control J2, T(| 11]. 



However, several problems related to the idea of modeling a flow by a small number of 
variables are open. One of the issues is the asymptotic stability of the models obtained. 
Often such models are capable of correctly reproducing the dynamics over small time 



intervals, whereas the asymptotic behavior converges to incorrect limit cycles II 1711 . 
This issue is related to both numerical artifacts and to an improper representation of 
the solution lfl5i Eol [Till . As a results low-order models are of delicate use and not 
robust to parameter variations. 

The present study describes a method to obtain robust low-order models. In previous 
works we showed that it is possible to obtain accurate low-order models of relatively 
complicated flows by minimizing the error between the model results and the reference 
solution J3]. Here we extend those works to cases where the flow is actuated by devices 
that can affect locally or globally the velocity and pressure fields. The objective is 
to derive a low-order model that provides accurate predictions and that is robust to 
variations of the control law employed. The main idea is to identify the manifold over 
which the non-linear dynamics of the POD modes lies, when the input to the system 
is varied. In this spirit, several dynamics are included in the identification procedure 
coupled with a Tikhonov type regularization. The case of a precomputed control as 
well as the case of a feed-back control are studied. 

The practical relevance of this work is that low-order models make possible to devise 
or to optimize controls for large-scale problems that would not be otherwise solvable in 
terms of computational size. Applications of this method is straight forward for models 
other than the Navier-Stokes equations. 



2 Reduced Order Modeling using POD 
2.1 Flow setup 

We consider a two-dimensional laminar flow past a confined square cylinder. This 
setup presents a reasonable compromise between physical complexity and computatio- 
nal cost. A sketch showing the geometry, the frame of reference and the adopted nota- 
tion is plotted in FigQ] At the inlet, the incoming flow is assumed to have a Poiseuille 
profile with maximum center-line velocity U. With reference to Fig[T] L/H = 1/8, 
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Fig. 1 - Computational domain ft. 




senses 



Fig. 2 - Placement of synthetic jet and sensors for control 



Li n /L = 12, L out /L = 20. No-slip conditions are enforced both on the cylinder and 
on the parallel walls. Details concerning the grids and the numerical set up are reported 
in 0|. 

All the quantities mentioned in the following have been made non-dimensional by L 
and U. The two-dimensional unforced flow obtained is a classic vortex street with 
a well defined shedding frequency. The interaction with the confining walls leads to 
some peculiar features, like the fact that the vertical position of the span-wise vortices 
is opposite to the one in the classic von Karman street @|. 

The presence of an actuator is modeled by imposing a new boundary condition on a 
small surface T c of dil : 

u(x, t) ■ n{x) = c(t), x gT c 

For control purposes we place two actuators on the cylinder. They are driven in opposite 
phase, as shown in Fig. [2]: 

v(x, t) = c(t), x G r c 

The control law c(t) can be precomputed, or obtained using a proportional feedback 
law. For example, using measurements of the vertical velocity at points Xj in the cy- 
linder wake, we can define a proportional control law : 

3 = 1 

where 7V S denotes the number of sensors used. 
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2.2 Proper Orthogonal decomposition with the snapshot method 

Seeking a reduced order solution that can be written 

N T 

UR.{x,t) = a r (t)<b r (x) 

r=l 

where the spatial functions <t» r are obtained by POD has become a popular approach 
when dealing with large scale systems. A vast literature concerning the POD procedure 
exists ll-ill . we refer to these works for a more general review of the method. 

2.2.1 The POD basis 

In our case, a numerical simulation of the Navier-Stokes equations is performed 
over a time interval [0, T], and the velocity field is saved at Nt time instants ti € [0, T]. 
This yields a data set {u l (x) = u(x,t l )}i = i,,N t . The aim of the POD procedure 
is to find a low dimensional subspace of C = spanfa 1 , ■ ■ ■ , u *}, that gives the 
best approximation of C. We therefore seek an orthonormal set {3> r } r =i...iv r , where 
N r << Nt, and a set of coefficients a\ such that the reconstruction error : 



N. 

U" — 

r=l 



(i) 

L 2 (n) 



is minimal. 

Following Sirovich's idea [21] the POD modes are expressed as linear combinations of 
the snapshots : 

N t 

$ k (x) = ^2b^u i (x) forfc = l.JV r 

2=1 

The vectors -iv ^ f° un d to be the eigenvectors of the Nt x Nt correlation 

matrix R, Rij = (u l , it- 5 ), corresponding to the N r highest eigenvalues, while the al 
are equal to the scalar products (u l , $ r ). 

In the case of forced flow, the snapshots depend on the control law c(t) used. In this 
work we consider POD basis derived from numerical simulations obtained using seve- 
ral different control laws. Indeed, there are a number of other parameters that could be 
varied, but since our aim is to study the effect of a control law, we set ourselves in the 
following framework : 

- Time instants, Reynolds number, domain geometry, placement of the actuators 
will be the same for all the snapshots in the database. 

- The control law c(t) will be varied 

The data set used for the POD is therefore written : 

{U h (x) =}i=l..N t ,t=l..N c 

where N c denotes the number of control laws considered. If C — {ci, C2, ■ ■ ■ , cn c } is 
the set of control laws used to obtain the database, the ensuing POD basis is denoted 
<J?(C). In the first part of this work, C is reduced to a single element which we denote 
c(t). 
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2.2.2 Dealing with the boundary conditions 

In the non-controlled case, we lift the boundary conditions on the velocity fields by 
defining a new set of snapshots : 

w' l (x) = u l (x) — u(x) 

where u is some reference velocity field that satisfies the same boundary conditions as 
the snapshots. In the present configuration, it can be the steady unstable solution, or a 
time average of the snapshots u k . 

When an extra boundary condition is imposed on the cylinder for control purposes, the 
snapshots are chosen to be : 

w l (x) = u l (x) — u(x) - c(t l )u c {x) 

where u c (x) satisfies the following criteria : 

u c (x) = u(x) on r\r c , u c (x) = 1 on r c 

In practice we use the velocity field proposed in (H : 

u c (x) = — (u'(x) — u(x)) 

c* 

where ft' is obtained in the same way as u but applying a constant control equal to c* 
on T c . The low-dimensional solution is now written : 

N r 

u R (x, t) = u(x) + c(t)u c (x) + ^2 a k (t)& k (x) (2) 

k=l 

2.3 POD-Galerkin Reduced Order Model 

Galerkin projection of the incompressible Navier-Stokes equations onto the first 
N r POD modes yields a system of ordinary differential equations : 

a r (t) = A? + C% r a k (t)+B% sr a k (t)a s (t)+P r 

+ E°c(t) + F°c 2 {t) + G°c{t) + H% r a k {t)c{t) 

(3) 

a r (0) = a° r 
1 < r < N r 
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where : 



((^•v)$ s ,* r ) 

Ct = -((u • V)# fe , <F) - ((* fe • V)u, * r ) + ^(A* fe , * r ) 
J3? = (u c ,# r ) 

G° = -((u • V)u c , *' r ) - (K • V)u, * r ) + ^-(Au c> *'•) 

F r G = (K-VK,$ r ) 

^° fe - (K • v)* fc , $ r ) + ((* fe • vx, * r ) 

7> r = (Vp,* r ) 

We note that since the snapshots satisfy the continuity equation, the modes do also. This 
implies that the pressure term V r is equal to j dn p$> r ds. If velocity field is constant at 
the boundaries, the POD modes are zero there. The pressure term therefore disappears 
completely. 
Setting : 



X r — A r , {B ksr } k s=1 ,„ Nr , {Cfer}fc=i...j\r r ! E r , F r , G r , {H kr } k=1 , 
and 

f(a(t),c{t),c{t)) = 1, {a k (t)a s (t)} k 8=v „ Nr , {a k (t)} k=1 ... Nr , 

c(t), c 2 {t), c(t), {a fe (*)c(*)} fe =i...ivj 
the first equation in (0 can be written in the compact form : 
a r (t) = f(a(t) 7 c(t),c(t))-X? 

The initial value problem (01 is a reduced order model of the Navier-Stokes equations, 
called the POD-Galerkin model. Such a model might be inaccurate for it may not take 
into account enough of the dynamics. Indeed, although a number N r of modes can be 
sufficient to capture most of the flow energy, the neglected modes continue to play an 
important role in the flow dynamics through their interaction with the conserved ones. 
The difference between the solutions a r (t ) of (|3j, and the coefficients a r (t) obtained by 
projecting the numerical data onto the POD modes (a r (t l ) = a*), has been underlined 
in several papers il 0] ■ It is therefore interesting to build a model that exploits the 
knowledge one has of the dynamics, that is the set of temporal projection coefficients 
a\ . This is the subject of the next section. 



3 Robust low order models 
3.1 Calibration method 

The idea of calibration is to keep the structure of the above model while adjusting 
the coefficients of the system so its solution is closer to the desired one. In previous 
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work 191 |8J], it was shown that robust low order models could be obtained by solving 
the minimization problem : 



min > / (a r (t) — a r (t)) dt 
subject to a r (t) = f (a(t) , c(t) , c(t)) ■ X r 



This state calibration method, which involves solving a strongly non-linear system, 
works well as long as the number of snapshots considered remains limited. For a large 
number of snapshots, the computational costs are excessive. 

Another method was suggested in and has been experimented, for a case with no 
control, in |0, 0] with good results. It consists in choosing X as the solution of : 

N r .T 2 

min^ / (h r (t) - f(a(t),c(t),c(t)) ■ X r ) dt (5) 

r=l ^° 

This method can be interpreted as approximating the error 

e r (t) = a r (t) - /(&(*), c(t), c(t)) • X G r 

by a quadratic function of all the non-discarded temporal coefficients, c(t) and c(t). 
Other choices for the approximation of e r lead to partial calibration problems. 
For example, if we suppose e r « Af. + Cf r ctk + G E c then we will solve : 



E / (firW - ' Ki - /aW • X°2) 2 dt (6) 

„ 1 J 



where 



mm 

"*1 r=l 



r,l — 



{Bksr}k,s=l—N T ' E r , F r , {-H& r }fc = i...j\r 



and 



/ 2 W 



!, {afc(*)}fe=i...jv r i C W] 

{a fc (i)a s (i)} fc V .. N , c(t), c{t), {a k (t)c(t)} k=1 . 



N, 



Of course, other choices of which terms to calibrate or not can be made. For a ge- 
neral formulation we denote N ca i the number of terms of vector X r that are calibrated, 
and we have iV ca ; ^ N% + 2 x N r + A. Whatever the choice for N ca i, this approach 
is always much more efficient than since it involves solving N r linear symmetric 
systems of size N^ al : 

ft{t)h(t) dt x c rS = fl(t) (k r (t) - f 2 (t) ■ x r G 2 ) dt (?) 
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The more terms of X r are calibrated, the more the problem becomes ill-conditioned. 
For this reason we choose not to calibrate the terms Bksr- 

Once the model has been calibrated to fit a particular control law c(t), it can of 
course be integrated using another control law. Denoting the input control law c'(t), 
the calibrated model is written : 



TZ({c}) 



a r (t) = f(a(t),d(t),ff{t))-X r 

Or(0) - a° r (8) 

1 < r < N r 



where by denoting lZ({c}) the model we put in evidence that it was calibrated using 
the control c(t). 

3.2 Well-posedness and robustness 
3.2.1 Calibration with feedback control laws 

We suppose that the control is obtained using a proportional feedback law (Sec l2.U : 

iV a 
3 = 1 

We can now consider two different calibration problems. The first is the problem (O, 
the second is : 



min^ f (h r (t)- f{a{t),c{t),'c(t))-X^ 2 dt (9) 

r=l •'° 

where c is defined by : 

c(i) = KjVRixj^) = Kj (^v( Xj ) + c(t)v c { Xj ) + arit^KxjU (10) 

This last approach makes the reduced order model a feedback model, which is useful 
if we want to use the model to determine an optimal feedback law. The problem is 
however under-determined. 

We reformulate dTOb to clearly show the dependency of c on a : 

N r 

c(t) = k + K r a r (t) (11) 
r=i 

where 

K ° = i - * {xj) ' Kr= i- w ^ 
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We now look at the partial-calibration problem described above. The function f 1 that 
appears in system (O can be reformulated : 

/*(*) = [!) { & k(t)} k=1 ... Nr , Ko + K t ae(t)] 

System (|7]i is therefore rank deficient. The problem remains if more of the system 
coefficients are calibrated, and according to the choice made the rank of the problem 
matrix can even diminish with respect to N ca i. 

This difficulty can however be solved by using one of the two methods proposed in the 
section ( 13.31 ). Finally, the proportional feedback reduced order model is written : 



d r (t) = f(a(t),c F (t),j(t))-X r 

C F (t) = £ K A + cF (*H(^) + E a r(tWv(*j 

j=l \ r=l 

a r (Q) = a° r , c F (0) = c(0) 
1< r < N r 



(12) 



3.2.2 Instability issues 

The system solved for calibration can be ill-posed even in cases different to the one 
just described. To understand why this is, it is sufficient to go back to the state cali- 
bration method mentioned at the beginning of the section. Solving the minimization 
problem (01 involves solving a non-linear system for which the uniqueness of solution 
is not guaranteed. The state calibration functional can therefore have several local op- 
tima, and so there are several possible choices for X that will lead to a low value of the 
error \\a — a\\. Since these choices should also be good choices for the minimization 
problem ©, the matrix L f[fi dt in (O is in general almost singular. A model obtai- 
ned by inverting this matrix is most often very unstable. To overcome this problem we 
propose a Tikhonov type regularization method which we describe in the next section. 



3.2.3 Robustness 

While a calibrated reduced order model lZ({c}) works well when integrated with 
c)(t) = c(t), its behavior when integrated with a different control law is unpredictable. 
As such, the reduced order model can difficultly be used for estimation and optimiza- 
tion purposes. 

In the literature several methods are proposed for adapting reduced order modeling for 
control purposes, some successful examples can be found in 

Odd. 

However for 

those cases, no calibration seems necessary for the models to work, but this is not the 
case for general control problems as shown in the following. 

The originality of the model we propose hereafter, is the combination of multi control 
data sets with the calibration procedure. Such a model is fast to build and yet remains 
accurate for different control inputs. 
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3.3 Building a robust low order model 

In this section we describe a method to make the reduced order model stable and 
robust. 

3.3.1 Tikhonov regularization 

Correctly solving (|4]i can be done by applying a quasi-Newton method, initiali- 
zed with a and Xf . It therefore seems reasonable to solve the following regularized 
problem, instead of (O : 

N r x 

where a is the regularization parameter. 

The parameter a can be chosen by a classical technique. We start by plotting, for a set 

2 

of values of a in [10~ 6 , 10~ 2 ], the error J2 r \\®r ~ *MI versus the coefficient varia- 
tion \\X° — Xf\\ 2 . This leads to a classical Tikhonov L-shaped curve of which the 
corner point is optimal in the sense that it is a good compromise between the error 



on the dynamics and the distance from the original coefficients [12]. The value of a 



corresponding to this point can be chosen to perform the calibration procedure. 



3.3.2 Calibrating over more than one control law 

In this paragraph we look at the changes to be made to the reduced order model 
when the data set includes simulations obtained using different control laws. Letting : 

the calibration problem becomes : 

N T N c T 2 
n ™EE / (^r(t)-f(& e (t),c e (t),c e (t))-X r ) dt (14) 
r=l e=i ^° 

We remark that although the size of the snapshots database is proportional to the num- 
ber of controls considered, the size of the calibration problem remains constant. Fur- 
thermore, if N c > 1 the rank deficiency discussed for proportional feedback no longer 
occurs. 

The main idea is that as the number of controls N c is increased, although the model can 
become a little less precise for the reference control, it is much more accurate for other 
control laws. In the next section we show some successful examples of this method at 
different Reynolds number, and for different kinds of control laws. 
We refer to a model built using 7V C control laws as an iV c -control model. Such a model 
is denoted TZc where C = {c\, ■ ■ ■ , cn c }. 
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4 Results and discussion 

The described technique was applied in order to build a low order model of the 
actuated flow around the confined square cylinder in various configurations. We tested 
the prediction capabilities of the model for two different Reynolds number, Re = 60 
and Re = 150, with precomputed and feedback control laws. In particular we built 
different models with one and more control laws and we analyzed their predictions 
with different controls. 

In all the examples presented in the following, actuation is started only once the flow 
is fully developed. With the control turned on the simulation is performed for about 
seven vortex shedding cycles, and N t « 200 snapshots are saved. T ~ 50 is the non- 
dimensional duration of the time interval. The number of POD modes retained for the 
reduced order model is N r = 40 for the case Re = 60 and N r = 60 for the case 
Re = 150. 

We measure the accuracy of the model 11(C) in the following way : 

- Time coefficients dynamics : 

For a given value of r, plot a r (t), solution ofH(C) with input c l (t), against a r (t), 
projection of the full order solution onto the POD basis 3>(C). In the examples 
r = 3 is usually chosen because it was the mode for which the differences bet- 
ween models were the most remarkable. 

- Computation of the integration error : 



In the examples with feedback laws we use only one sensor placed in the cylinder wake. 
Choosing the center of the cylinder as the origin of a coordinate system, we denote 
x s = (x s ,y s ) the position of the sensor. The integration error £ F (C, K 1 ) is measured in 
the same way as for the non-feedback case. 

Our first goal is that the model should be able to reproduce the DNS data to which 
is was fitted, we therefore expect £(C, c 1 ) to be small if c 1 G C. Our second goal is 
that the model be robust to parameter variation. As the difference between c'(t) and 
the controls in C increases, the error £(C, c') grows. We seek a model for which this 
growth rate is as low as possible. 

4.1 Divergence of a classical Reduced Order Model 

A simulation at Re = 60 was performed using feedback control with a sensor pla- 
ced at (x s , ys) = (0.7, 0.0) and K = 1. We denote c(t) the control law obtained at the 
end of simulation. 

We compare the results obtained with the POD Galerkin model <[3j and with the cali- 
brated model 1Z({c}) (see system (0 for model formulation). The model integration 
error £({c}, c) is equal to 23% for the non-calibrated model, and to 0.136% for the 
calibrated model. 

For a feedback model, the difference is even more important. We integrated the feed- 
back system ( fT2b with K = 1, once with X obtained by Galerkin projection, and once 



£(C,c l ) = \\a-a\\/\\a\\ 
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FlG. 3 - Projection of the DNS simulations onto POD modes vs. integration of the 
dynamical system dT2T > with X = X G 
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Non-dimensional time Noi Mlhiensional time 



FlG. 4 - ci (continuous line) versus c F , when the model is calibrated (left) and when it 
is not (right) 

with X calibrated as described in 13. 2. II We obtained an integration error £ b ({c}, 1) of 
117% in the first case, against an error of 4% in the other. An example of the errors in 
terms of time dynamics that the non-calibrated model can produce are shown in Fig[3] 
In Fig|4]we plot the control law c F (i) computed when integrating the feedback model, 
and on the same figure, the original control law c(t). Results for the non-calibrated case 
are plotted on the right : the distance between c F (t) and c(t) increases with time, mea- 
ning that at each time step, new errors are added to the model. Calibration is therefore 
all the more essential when considering feedback control. 

In order to calibrate, regularization is needed to get well -conditioned inverse problems 
as shown in the following. However, the choice of the parameter a is not an easy one. 
For example, we performed a simulation at Re = 150 using a feedback control with 
a sensor placed at (x s ,y s ) = (0.7,0.0) and K = 0.8. The calibration described in 
|3.3.1| was performed with a ~ 0. This led to an ill-conditioned system to solve and to 
a model that was not very accurate, and not robust at all to parameter variations. The 
effect of a on model results is shown in Fig |5] The two top figures show the third mo- 
dal coefficient obtained by projection and by integrating the model with K = 0.8. On 
the left, we plot the results obtained when the model was built with a = 1.6 * 10~ 6 : 
at the end of the time period the model diverges from the DNS results. With a higher 
value, a = 10~ 3 , this problem no longer occurs, as shown on the right. The same test 
was then performed with a different value of K in order to see the models capacity to 
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Fig. 5 - an DNS (continuous line) with K = 0.8 (top) and K = 1.3 (bottom) versus 03 
obtained when the model is calibrated with a = 1.6 * 10~ 6 (left) and when a = 10~ 3 
(right) 

predict dynamics to which it was not fitted. The results are shown in the same figure : 
divergence was immediate for a low value of a, whereas for a higher value, the mo- 
del, although not accurate, was at least stable. It appears that, when using an higher 
regularization parameter, the calibration system is well conditioned, the model more 
accurate and more stable when integrated with a different control law to those used for 
calibration. In the following the parameter a is determined using the L-method with 
the restriction that any values of a below a certain threshold are excluded. 

4.2 Testing model robustness 

In this section we present the improvements brought to model robustness by intro- 
ducing calibration over several control laws. For both Reynolds numbers, Re = 60 and 
Re = 150, the same experiment was performed : 

Step 1 : Build 1- , 2- and 3-control models 

We started by choosing three control laws which we denote Ci(t), C2(t) and c${t). 
For each control we performed a simulation of the Navier-Stokes equations, saving 200 
snapshots for each simulation. We then defined seven control sets : 

Three 1 -control sets : C 1 = {ci},C 2 = {c 2 },C 3 = {c 3 } 
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Three 2-control sets : C 4 = {ci, C2}, C 5 — {01,03}, C 6 = {02,03} 

One 1 -control set : C 7 = {c\, c 2 , 03} 

For each control set C\ we computed a POD basis &(C l ) as described in section [2". 2. II 

and a calibrated reduced order model 1Z(C 1 ) by solving problem (fl4l) . 

In the following we refer to c\ (i), 02 (t) and 03 (t) as the model control laws. 

Step 2 : Run the model with different control laws 

We next chose several other control laws which we denote c* es '(£). Each of these 
test control laws was used as input for the Navier-Stokes equations, and for the seven 
reduced order models 1Z(C 1 ) described above. The snapshots from the Navier-Stokes 
simulations were projected onto the seven POD bases 3>(C l ). This procedure made it 
possible to compute the model integration errors £| = £(C l , c* est ), and compare the 
efficiency of each model. 

For measuring model robustness, it is useful to have some idea of how much the dy- 
namics we are trying to predict, differ from those included in the model. We therefore 
need to find a way, for each model, to measure the distance between the N t x N c snap- 
shots that were used to build it, and the N snapshots obtained using a test control law. 
To do this we proceed in the following way : if the control set C l is composed of N c 
control laws, then the distance between the simulations associated to C\ and the one 
obtained using c* e;st (i), is defined as : 



where the terms a n (n~jorn— 1 ■ ■ - l) result from projecting the snapshots onto the 
POD basis $(C J ). 

The results are plotted for in Fig. [8] and Fig. Q~2] For each value of model i, the model 
integration error E\ is plotted versus the distance A? . We note that the three controls 
used to build the models were in fact included in the test set, which explains why there 
are 3 points at = 0. 



4.3 Results for Re = 60 

In Fig|6]we plot the control laws used to build the models. For each control law we 
plot the third modal coefficient 03 (t) to give an idea of the dynamics induced. The fi- 
gure also shows the prediction for this coefficient given by the 3-control model 1Z(C 7 ). 
The model results are accurate : the reduced order model was successfully calibrated to 
fit several dynamics. Eleven extra control laws were used for testing. A few examples 
are plotted in Fig. [7] For these examples we also plot the third modal coefficient obtai- 
ned by projection and by model integration. Some discrepancies in coefficient ampli- 
tude are observed, but overall the model predicts the right time dynamics. 
In Fig.|8]we look at the results obtained with the different models, using the distances 
and errors described above. The first point to be made is that the model error is almost 
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FlG. 6 - Control laws used to build the models (top) ; 03 DNS (continuous line) vs 
prediction by 3-control model (symbols) 



zero when the distance from the model is zero. This confirms that 1 -control models 
work well when integrated with the control law to which they were fitted. The errors 
then increase with the distance from the model, as was expected. 
The graph highlights the disadvantage of 1 -control models. In the best case the diffe- 
rence between projection and prediction coefficients becomes higher than 20% as soon 
as the distance from the model exceeds 40%. In contrast, for the 2-control and 3-control 
models, the error stays under 20%, even when the distance increases. In Fig. [9] we plot 
isolines of the vorticity at time t = T for one of the test control laws (the third control 
law in Fig[7j. Time coefficients were obtained by solving 1Z(C) with C = {ci,C2}. 
The velocity field was then reconstructed using the first ten of these coefficients and 
the first ten POD modes in 3>(C). The reconstructed vorticity is presented along with 
the vorticity obtained by running the Navier-Stokes equations with the test control law. 
The controls used to build the model caused a slight decrease in vortex size (see Fig. 
[6] bottom left) whereas actuation used in the test caused a slight increase in vortex size 
(see Fig. [7] bottom right) . We note that the model was able to predict such features, 
and that at the end of the simulation time, the structure of reconstructed flow is almost 
identical to that of the real flow. In contrast, the 1 -controls weren't able to identify this. 
If the same reconstruction is performed using C = {ci} for example, the flow appears 
almost stable at t — T, meaning the model predicted the opposite behavior to what 
actually happened. 
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Fig. 7 - Control laws and time coefficients used for testing the model 



80 




Distance (%) 



FlG. 8 - Prediction errors obtained using 1 -control, 2-control and 3-control models 
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FlG. 9 - Model predicted vorticity field (top) and Navier-Stokes vorticity field at t = T. 
Positive (continuous lines) and negative (dashed lines) vorticity isolines 

4.4 Results for Re = 150 

For Re = 150 only feedback control laws are used both to build the models and to 
perform the tests. In Fig[lO]the three feedback control laws used to calibrate the model 
are shown. The laws are obtained with sensors placed at (x s ,y s ) = (0.7, 0.0) and by 
using gains K = 0.6, K = 0.8 and K = 1. The figure also shows the third modal coeffi- 
cients given by integrating the 3-control feedback model with each gain. Although the 
control laws induce different dynamics, the model is able to give an accurate predic- 
tion in all three cases. Six extra control laws were used for testing, each corresponding 
to a different choice of K. A few examples, with corresponding coefficients 0,3 (t) are 
plotted in FigQT| It appears that the dynamics are quite different when the distance, 
between the gain value and gains included in the model, is large. For example, when 
using a gain K = 0.1 the average value of 03 (i) is low compared to that obtained 
with K = 1. However, the 3-control model again gives an overall good prediction of 
the time dynamics. Fig[T2] is built the same way as Fig. [8] In particular the graph 
shows the disadvantage of using a 1 -control model, with prediction errors of over 34% 
when the distance from the calibration dynamics increases over the 30%. As in the 
case Re = 60, the 2-control models give more accurate predictions than the 1 -control 
models. The lowest errors are obtained with the 2-control model (K = 0.6, K = 0.1). 
This observation suggests that, in model construction, an optimized a priori choice of 
the sampling points could be useful to obtain a more robust model. We note that in this 
case it was the model built to fit the highest and lowest values of K that gave the best 
result, and that adding a third intermediate control to the model (K = 0.8) did not bring 
any improvement : the 3-control model gives more or less the same results. 
In Fig. Qj] we plot isolines of the vorticity at time t = T for the flow obtained using 
K = 0.1 as feedback gain (the first one in FigfTTh. Time coefficients were obtained by 
integrating the 3-control model. The velocity field was then reconstructed using all the 
60 coefficients and POD modes. The reconstructed vorticity is presented along with 
the vorticity obtained by running the Navier-Stokes equations with the test control law. 
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FlG. 10 - Control laws used to build the models (top) ; 03 DNS (continuous line) vs 
prediction by 3-control model (symbols) 

The controls used to build the model were similar in the sense that they had a much 
stronger effect on the flow compared to the control obtained with K = 0.1. We note 
that the model is able to accurately predict a flow snapshot and that the reconstructed 
flow is almost identical to that of the real flow. 

5 Conclusions 

The overall picture of reduced-order modeling that results from our study is the fol- 
lowing. Given a control law, one can deduce a low-order model of the actuated flow by 
simply projecting the Navier-Stokes equations on POD modes. The coefficients of the 
quadratic model thus obtained are found by projection. However, a model constructed 
this way will show large time-integration errors even for the same control law used to 
generate the POD modes. Calibration can take care of that, in the sense that the mo- 
del coefficients can be determined in order to match as closely as possible at least the 
solution from which the POD modes are obtained. This might lead to a numerically 
stable model. However, this model is generally not at all robust, in the sense that the 
predictions for a slightly different configuration from that it was generated from, fails. 
A symptom of such lack of robustness is observed in the ill-posedness of the inverse 
problem : the matrices to be inverted are almost singular. 

In order to get around this deficiency, we regularize the solution by adding a constraint 
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Fig. 11 - Control laws used to test the models (top); 03 DNS (continuous line) vs 
prediction by 3-control model (symbols) 

to the minimization method used to solve the inverse problem. We ask that the coef- 
ficients of the polynomial expansion be close enough to those obtained by projection. 
This method allows to synthesize models that adequately simulate the flow in a small 
vicinity of the control law used to generate the solution database. However, the actual 
real improvement in robustness is obtained by spanning the solution manifold, i.e., by 
including several control laws in the inverse problem definition. By doing this, the re- 
sults presented show that the models are able to predict dynamical behaviors that are 
far, in terms of an energy norm, from the cases included in the database. A consequence 
of such an additional regularization is that the matrices involved in the inverse problem 
solution become well conditioned. 

Another important aspect of the method proposed, is that its cost is that of a matrix in- 
version, and that it does not scale with the number or the size of data sets used to build 
the model. Therefore it seems reasonable to envisage an automatic strategy to enrich 
the model by spanning the control space. In this respect, the technique proposed in |5] 
to distribute in an optimal way the points where to test the control space can help mi- 
nimize the number of a priori simulations needed to build the model. For example, our 
results show that a model based on two controls might predict the effect of actuation 
laws not present in the data base, as precisely as a model based on three controls, if the 
two controls are appropriately placed. 

In conclusion, the modeling we propose appears to be a viable approach to determine 
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FlG. 12 - Prediction errors obtained using 1-control, 2-control and 3-control models 




Fig. 13 - Model predicted vorticity field (top) and Navier-Stokes vorticity field at 
t = T obtained with K = 0.1. Positive (continuous lines) and negative (dashed lines) 
vorticity isolines 

control strategies for those problems that because of their computational size cannot be 
treated in the framework of classical control theory. 
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